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Single impurities in a Bose-Einstein condensate can make two polaron flavors 
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Polarons, self-localized composite objects formed by the interaction of a single impurity particle 
with a host medium, are a paradigm of strong interaction many-body physics. We show that dilute 
gas Bose-Einstein condensates (EEC's) are the first medium known to self-localize the same impurity 
particles both in a Landau-Pekar polaron state akin to that of self-localized electrons in a dielectric 
lattice, and in a bubble state akin to that of electron bubbles in helium. We also show that the 
BEC-impurity system is fully characterized by just two dimensionless coupling constants, and that 
it can be adiabatically steered from the Landau-Pekar regime to the bubble regime in a smooth 
crossover trajectory. 

PACS numbers: 05.30.Jp, 03.75.Hh, 67.85.Hj, 67.85.Bc 



The polaron, a single distinguishable particle that 
interacts with the self-consistent deformation of the 
medium that contains it, is a paradigm of strong in- 
teraction physics in condensed matter [1] [2], chemistry 
[3], and biophysics [4j. Polarons self- localize when suffi- 
ciently cold and strongly coupled to the host medium. In 
nature, large [5^ polarons appear in two flavors: particles 
that hardly deform the medium, such as electrons in di- 
electric lattices |61 , and particles that greatly distort the 
medium, such as electron bubbles in condensed helium 
superfluids [7]. In the first class, the particle is accompa- 
nied by a cloud of small amplitude collective excitations 
of the medium which, as first shown by Landau and Pekar 
[Sill], can be integrated out in the strong couphng hmit. 
In the bubble systems, which occur in fiuids and dense 
gases, Kuper [10] showed that the strongly repelling par- 
ticle can be described as residing in a self-created cavity 
- the "bubble" - surrounded by the fiuid. The effective 
mass and mobility of bubble and Landau-Pekar polarons 
exhibit quite different behaviors, and so the two polaron 
flavors are customarily treated as distinct. However, 
as we show below, dilute gas Bose-Einstein condensates 
(BECs) [m [12], unlike previously known host media, 
can self-localize neutral impurity atoms [131 [H] in both 
the bubble and Landau-Pekar polaron states. Further, 
the system phase diagram presented here shows how the 
strongly coupled BEC-impurity polaron evolves contin- 
uously between these limits as the interaction strengths 
and BEC density are varied. The BEC-impurity system 
can therefore be regarded as a quantum simulator of large 
polarons in a boson environment I5j, complementing re- 
cent proposals for simulating lattice polarons [T5H19) . 

The system: We consider a neutral impurity atom of 
mass TO/ immersed in a homogeneous BEC of N boson 
particles of mass ms contained in a macroscopic vol- 
ume J7, giving an average density p = N/fl. Bosons at 
positions r and r' interact via a repulsive short-range 
interaction of scattering length aBSi described by an 
effective potential Vbs (r — r') = XbbS {r — r'), where 
Abb = {4:TrH'^ /ms) ass and ubb > 0. The boson chemi- 
cal potential /is ~ ^bbP sets the time scale h/ ps and the 



coherence length ^ = l/^/lGirpaBB sets the length scale 
on which the BEC can respond to a perturbation. The 
perturbation here is provided by the impurity interacting 
with the bosons via Vjb (r — x) = XjbS (r — x), where x 
represents the impurity position and the impurity-boson 
interaction strength Xjb ~ 2TTh'^ (mj + m^ ) ajB is 
proportional to the impurity-boson scattering length uib , 
taken to be Feshbach-tuned ^20] to a large positive value 
[H] . We break the translational symmetry of the BEC- 
impurity system ground state by fixing the impurity 
center-of-mass position at r = [12]. We write the 
BEC density in the presence of the impurity as Pb (r) = 
p + 6pB (r)- Self-localization occurs when the effective 
potential Xib^pb (r) can trap the impurity. We classify 
\SpB{y^ = 0)/p\ < 0.1 as the Landau-Pekar regime [T^ . 
and \Spb{tc = 0)/p\ > 0.9 as the bubble regime. We will 
show that in these regimes the impurity observables ex- 
hibit the scaling behaviors expected from the simplified 
Landau-Pekar and bubble descriptions respectively. 

Landau-Pekar polaron: The impurity-boson repulsion 
can be simultaneously strong enough to self-trap the im- 
purity, and weak enough to hardly change the BEC den- 
sity profile [Fig.jlja)]. In this regime, the Bogoliubov ex- 
pansion and transformation that describes the BEC fluc- 
tuations |23| can be carried out around the homogeneous 
BEC [M]- The interaction of the dilute BEC with an 
impurity of density p/ (r) = (27r)~ J d'^ke'^^'^pi^k then 
gives a Frohlich hamiltonian [25], familiar from electron- 
phonon interactions. Representing the quasi-particle an- 
nihilation/creation operators of momentum k and energy 

hiuk = hk^/iiTeWil^J^ by 

boson interaction is described by 



^k,^k) the impurity- 
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^Pl~k PBM 
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A/BP + -/^Yl ^kPi,-k [bl + 6kj , (1) 



where /OB,k is the operator associated with the boson den- 
sity. In Eq. (|T]), the impurity-phonon interaction matrix 
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FIG. 1; Numerical results for the normalized boson density 
(blue shading) and un-normalized impurity wavefunction (red 
dashed line). Parameter values [see Eqs. (|3| and Q]: (a) 
/3 = 25 and a = 10-^ (b) /? = 25 and a = 10-^ (c) /3 = 25 
and a = 10^ (d) /3 = 5 x 10'' and q = 10^ 



element M = A/i3^ , and Vk = {^^P/ [l+^^k^])^^* 
is a structure factor arising from the Bogoliubov trans- 
formation. By showing that the Frohlich coupling Hjg 
amounts to a displacement of the k-mode oscillator co- 
ordinate 



^i. - {bl + 6_k) /\/2, 



Landau and Pekar |9i 



integrated out the phonon modes. The resulting energy 
reduction takes the form of a self- interaction [531 IlS]j 
AE = -M2(27r)"^/d3kpj_kP/,k(i'fc/Swfe), which in 
the strongly coupled regime overcomes the kinetic energy 
cost of localizing the impurity. Here, the self-interaction 
potential is an attractive Yukawa potential of range £,. It 
can cause self-localization when ^ exceeds the impurity 
extent, which is comparable to |13j 

Ro = ['iTTpajg (1 -I- mj/mB) (1 -l- mBlmi)] . (2) 

Specifically, when the ratio of the boson healing length ^ 
to the self-localization length R^ 
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the above description predicts self-localization jTB] . Note 
that the impurity extent Ro — ^/ P can be significantly 
smaller than ^, making this polaron a sub-coherence 
length structure. The binding energy, proportional to 
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mj 
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(4) 



then significantly exceeds fiB- When scaled by Eg and 
Ro^ impurity observables in the Landau-Pekar regime de- 
pend only on the dimensionless coupling strength /3 |13j . 



Bubble polaron: When ajB grows sufficiently to ex- 
pel the EEC from the impurity's vicinity [Fig.nTc) and 
(d)], a Bogoliubov procedure should expand around the 
deformed EEC. However, a bubble description such as 
Kuper's model of electron bubbles in helium is a sim- 
pler starting point 10 . With complete BEC and im- 
purity separation [Fig.flTd)] the impurity, trapped in a 
self-created spherical cavity of radius Re and volume Vc, 
has wavefunction x (^) = \/t^Rc sin (jxr j R^) jr. Ne- 
glecting surface tension, the system energy difference E^ 
with and without impurity is the impurity kinetic energy 
-K^T? I (mjR^^ plus the energy cost of making the cavity, 
PVc, where P ~ XbbP^ /2 is the BEC pressure. Hence 
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(5) 



mjR'^ 3 TUB 

The minimization dEc/dRc = yields the expected 

cavity radius Re = [4:{mj/mB) P^ubb] , and 
the impurity energy Ec = (5/3) (tt^S^/ [777/i?^]) = 

(57r/2iV5) (m^/mjf^pB/ [y/^^f^'- 

BEC permeability: The radical change in BEC- 
impurity overlap seen in Fig.llTa) - (d) is due to a BEC 
"stiiTness" arising because it costs kinetic and interaction 
energy to move ANb = |/rf'^r [pb (r) — p]\ bosons away 
from the impurity. Estimating the interaction energy 
cost as E^ = ANbPb and using the predicted displaced 
boson number from [57], ANb = \^ib/^bb\ (valid in 
the Landau-Pekar and crossover regimes but not in the 
bubble regime), gives E^ — \\ib/^bb\pb- The ratio 

mi Y I tub ' ^ 
rriB J \mj 
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then quantifies the relative importance of the displace- 
ment energy cost to the overlap energy gain of self- 
localization. We refer to a as the BEC-impurity "perme- 
ability" . In the Landau-Pekar regime, a direct analytical 
evaluation yields \SpB (r = 0) /p| — (4-\/2/3v^) <t^^ = 
1.064 cr~^. Thus, tJ ^ 1 implies Landau-Pekar conditions 
where the repulsion is insufficient to overcome the BEC 
stiffness and displace the bosons noticeably [Fig.fTJa)]. 
A gradual increase in a/s then expels the bosons sig- 
nificantly when a ^ 1 [FigJIJb)], and enters the large- 
depletion bubble limit [Fig.fTTc) and (d)] when cr ^ 1. 

General case: A more general ground state treatment, 
encompassing the Landau-Pekar and bubble regimes as 
limits, is based on the strong coupling approximation of a 
many-body product state. Minimizing the energy while 
requiring the respective boson and impurity wavefunc- 
tions ^ (r) and x (r) to be normalized gives two coupled 
Gross-Pitaevskii equations 

^'^'v(r) + ABB|V'(r)N(r) 



fJ-Bip (r) 



e/X (r) 



2777s 

+A,s|x(r)N(r) 



(7) 



2777/ 



x(r) + A/B|V(r)rx(r). 



where limr-j-oo '0(^) = \fp^ ^ind \ib and e/ repre- 
sent the Lagrange multiphers ensuring normahzation, 
/d^r |V.(r)|^=7V,/d3r |;t(r)|2 = l. 

The BEC-impurity system has five physical parame- 
ters - mB,rni, p,aBB and o/b - but we find that the 
proper dimensional scaling of energies, density, and dis- 
tances reveals a minimal dependence on just two coupling 
constants. The first is the length scale ratio /? = ^/Ro- 
The second is the mass-scaled boson gas parameter 
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P^BB 
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All properly scaled observables can be cast in terms of 
a and /3. For example, the permeability parameter a 
takes the form a {a, 13) — l/(47r^'^-\/a^), and in the 
bubble limit the energy and cavity radius are given by 
E^ = Eo(J>TT/2A)p'^{A/af'^ and Re = Ro'^i^^^'^pa^'^ , 
respectively. To prove the minimal dependence, we sub- 
stitute (real-valued) scaled dimensionless boson {p) and 
impurity {g) wavef unctions; 



■0 (r) = 7^ p (x = r/i) 
X(r) = i?-3/2g(y = r/i?„) , 



(9) 



into the coupled equations ([7]), scaling the energies and 
lengths to obtain the dimensionless form 



l-p2(x) 



p(x) 



47r/32 
a{a,[i) 



g2 (y = x/3) p (x) 



( V", + ei) g (y) ^ ^a [a, /3) p^ (x = y//3) g (y) , (10) 

where ej = ej/Eo, and the g,p solutions have to satisfy 
the normalization J (Pyg^ (y) — 1 and boundary condi- 
tion lima:_>oo P (x) = 1. Since the scaled coupled equa- 
tions depend only on a and /3 |28] the BEC-impurity 
phase diagram is two-dimensional (2D). 

Numerical solutions of the coupled equations {Jj pro- 
vide a rigorous test of the Landau- Pekar and bubble lim- 
iting behaviors by determining the a- and /3-dependences 
of the relevant observables. For example, in Fig. [2] 
we plot the root mean-square impurity extent Rrms = 

Jj (Prr'^pi{v). Each curve shows Rrms/Ro as a func- 
tion of a for a specific /3-value on a log-log plot, with 
arrows indicating a values corresponding to unit perme- 
ability a — \. On the left, the a-independence of the 
Landau-Pekar impurity properties is confirmed by the 
zero slope of the curves, which also show the expected 
convergence to Rrms ~ 4.6 i?o for /3 > 20 [TT. On the 
right, the straight lines of slope 1/5 confirm the expected 
bubble scaling Rrms/Ro oc j3a^/^. The smooth change 
between these limits indicates a crossover, rather than a 
transition, between Landau-Pekar and bubble regimes. 

Likewise, the relative BEG density change at the impu- 
rity position, \5pB (r = 0) /p|, plotted in Fig.jsjshows the 
CT~^ ex a^'^ scaling of the Landau-Pekar polaron on the 
left side. Before crossing the line of unit relative density 
response, the curves level off and approach the bubble 
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FIG. 2: Dependence on gas parameter a of the r.m.s. width 
of the impurity density (in units of scaled length _Ro), for 
several values of impurity-BEC interaction parameter /3. For 
each P arrows indicate the values of a for which permeability 
a = 1. 
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FIG. 3: Dependence on gas parameter a of the relative 
decrease in BEG density at the impurity location, for several 
values of impurity-BEC interaction parameter /?. For each /3 
arrows indicate the values of q for which permeability a = \. 



limit of maximal BEG depletion asymptotically. As in 
Fig. [2] the arrows on Fig. [3] indicating the u = 1 posi- 
tions illustrate that the crossovers occur near unit per- 
meability. Fig.|4] is the 2D phase diagram for the BEG- 
impurity system, produced by plotting \8pb{^ = 0)/p| in 
the (a, /3)-plane, colored to show the polaron regimes and 
the region where the impurity-BEG interaction is not 
strong enough to cause self-localization. The red a = \ 
line lies, as expected, on top of the crossover region. 

Experimental realization: The remarkable ability of 
a BEG to self-localize impurities both in Landau-Pekar 
and in bubble states could be strikingly illustrated by an 
experiment that adiabatically Feshbach-tunes the same 
BEG-impurity system from one limit to the other. For 
typical densities (10^^ cm~'^ < p < 10^'' cm~^) and re- 
alistic ranges for a^s and mB/rni, we find that cold 
atom a values may range from 10~^ to 10~^. Increasing 



50 



i 40 



.9 -Tt 
■P jO 



20 



10 



1 1 
o-=l 1 

1 




I 

Landau-Pekar i 


Bubble 


Polaron 


, Crossover Polaron 


. (<5pg/p<0.1) 


• (Spb/p > 0.9) 




• 

1 


: i 


» 

» 
* 


Delocalized 

1 1 1 1 



10-10 10-6 10-2 

mass-scaled gas parameter a 



102 



FIG. 4: Phase diagram of the BEC-impurity system obtained 
as a plot of |(5ps(r = 0)/p|. The dots in the Landau-Pekar, 
crossover, and bubble regions correspond to the BEG density 
profiles in Fig.ma), (b), and (c), respectively. The blue line is 
a trajectory that tunes the system continuously from Landau- 
Pekar limit to crossover region to bubble limit (off the top of 
the plot). 



aiB by Feshbach tuning could achieve self-localization at 
/3 ~ 5 (see Fig.[4|. Using Eq. ([s]) and scaling scattering 
lengths and densities by the typical values of 1 nm and 
p — VS^^ cni~^, this corresponds to a/s = af;^-, where 



s.l. 



168 nm 



I Pl-p 



1/4 



(11) 



This self-localization near /3 ~ 5 results in a Landau- 
Pekar polaron if cr = [47r^/*-\/l25a] ^ 1, requiring 



1 /Q 

osB <C 2.00nm(p/p) {mB/mi) 



.-2/3 



(12) 



A further increase in ajs and/or p can lower the perme- 
ability and effect a cross over to the bubble regime when 
Or™**, where from Eq. (l6| 
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aJ'B°"' = 200 nm 



(P/P) ^^^(mB/mi) ^/^ 
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Throughout, the permeability parameter varies as 
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where afg denotes the scattering length of Eq.(ll) at 
standard density, p — p. The adiabatic {aiB,p) varia- 
tion traces out a trajectory on the phase diagram (e.g. 
blue line on Fig.[4| that starts in the low a, cyan-colored, 
Landau-Pekar region. An increase in BEC density (cor- 
responding to an exponential increase in (In a, /?)- space) 



and a Feshbach increase in ajB (corresponding to a ver- 
tical upward motion in the same space) can then, even- 
tually, steer the BEC-impurity across the green-colored 
crossover regime into the yellow-colored bubble region. 

We now consider two potential issues facing such an ex- 
periment. The first is the lifetime of the impurity against 
three-body recombination. An increase in ajB is gen- 
erally (but not always - see below) accompanied by a 
decrease in lifetime ^\. Estimating the three-body lim- 
ited impurity lifetime r/ as in (3D], we expect rf'** ~ 

(V3/3.9) a/1 -h 2 {tub /mi) (h/mB) ajgp^ . As the 

time scale of the BEC response (and the slowest time 
scale in the system), we expect tb — ti/pB to set the 
scale of the self-localization dynamics. The estimated 
impurity lifetime with full overlap can then significantly 
exceed tb as long as /3 is not too large. 



TB 



47r2 V3 (1 + mi/ruBf (1 + mB/mif 



/32 3.9 



a/1 + 2mB/mj 



(15) 



A more careful study of the three-body loss |3T] found 
that rf^^ should be divided by an oscillating Stuckelberg 
factor related to three-body Efimov physics [32]. Near 
the nodes of the Stuckelberg factor the lifetime greatly in- 
creases. Further, the near complete separation of the im- 
purity and the BEC in the bubble limit implies a signifi- 
cant increase in impurity lifetime for sufficiently large ajB 
because the three-body recombination loss rate is propor- 
tional to the overlap / d^r p\ (r) pi (r). We note that in 
condensed helium, the increased lifetime of positronium 
[331 [31] has been used as a signal of self- localization [33] . 

The second challenge faced by an experimental real- 
ization of the bubble polaron may be the buoyancy force 
F = -VEc (p) = - (4/5) Ec (Vp) /p attempting to expel 
the impurity from the high density region in an inho- 
mogeneous BEC. A two-color trap or a species specific 
potential [361 137| may be necessary to keep the impurity 
bubble in place. Alternatively, a homogeneous conden- 
sate might be created in a flat-bottomed Painted Poten- 
tial [3B]. We note that while the bubble polaron always 
seeks low BEC density, the Landau-Pekar polaron can be 
high density seeking. 

Conclusion: We have shown that a distinguishable 
neutral atom embedded in a dilute BEC can, if the BEC- 
impurity repulsion is strong enough, self-localize in a 
bubble polaron state in which the impurity is imperme- 
able to the BEC host medium. This state is analogous 
to that of an electron bubble in condensed helium. Re- 
markably, and uniquely, the same BEC medium can also 
self-localize impurities in Landau-Pekar polarons if the 
BEC density and interaction strengths have appropriate 
values. We characterized the overlap of the self-localized 
impurity with the host fluid by a permeability parame- 
ter (T that ranges from cr 3> 1 in the Landau-Pekar limit 
to (T ^ 1 in the bubble limit. We have shown that the 
BEC-impurity system is fully characterized by just two 
dimensionless coupling constants. In the corresponding 



phase diagram, the bubble and Landau-Pekar states cor- 
respond to broad regions that are separated by a smooth 
crossover region near a ^ 1. Finally, we pointed out 
that a single impurity-BEC system can be adiabatically 
steered from one limit to the other. 
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